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'• Abstract 

We study 12 parameter families of two qubit density matrices, arising from a special class of 
two-fermion systems with four single particle states or alternatively from a four-qubit state with 
amplitudes arranged in an antisymmetric matrix. We calculate the Wooters concurrences and 
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^~^, the negativities in a closed form and study their behavior. We use these results to show that 
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distributed entanglement. An explicit formula for the residual tangle is also given. The geometry 
of such density matrices is elaborated in some detail. In particular an explicit form for the Bures 
metric is given. 



I. INTRODUCTION 

Entanglement is the basic resource of quantum information processing^. As such it has 
to be quantified and its structure characterized. For entanglement quantification one uses 
special classes of entanglement measures which are real- valued functions on the states. Pure 
and mixed state entanglement and its quantification in its bipartite form is a well-understood 
phenomenon. Moreover, on the geometry and structure of entangled states associated with 
such systems a large number of interesting results is available^. 

For example for pure states of bipartite systems the classification of different entanglement 
types is effected by the Schmidt decomposition. If the Schmidt decomposition is known, 
from the Schmidt numbers one can form the von-Neumann entropy^ as a good measure 
characterizing bipartite entanglement. For quantifying mixed state entanglement no such 
general method exists. For the special case of two qubits as a measure of entanglement we 
have the celebrated formula of Hill and Wootters- for the bipartite concurrence C and the 
associated entanglement of formation. The structure of this measure of entanglement was 
studied in many different papers^. Its structure has been related to antilinear operators^, 
combs and filters^, and has also been generalized to rebits^. Explicit expressions for different 
special classes of density matrices and a comparison with other measures of entanglement 
has been given^i^iii. 

In this paper we would like to study the structure of special 12 parameter families of 
two-qubit density matrices for which the the mixed state concurrences can be calculated in 
a closed form. Such density matrices can be regarded as reduced ones coming from some 
larger system with special properties. In order to motivate our choice for this larger system 
we consider an example. If we consider a three-qubit state {ip) G C^ (8> C^ (8> C^, then after 
calculating any of the reduced density matrices e.g. Qu = TT3{\il'){'ip\) we are left with a 
two-qubit density matrix of very special structure. For example in this case Qu is of rank 
two, and this observation enables an explicit calculation of the mixed state concurrence in 
terms of the amplitudes of the three-qubit state {ip)- This result forms the basis of further 
important developments namely the derivation of the Coffman-Kundu-Wootters relation of 
distributed entanglement^^. 

Proceeding by analogy we expect that four-qubit states of special structure might provide 
us with further interesting examples of that kind. Let us consider a four-qubit state |\I^) G 



C^ C^ ® C^ ® C^. A class of two-qubit density matrices arises after forming the reduced 

density matrices like Q12 = Tr34(|\l/)(\1'|). However, density matrices of that kind are still 

too general to have a characteristic structure. Hence as an extra constraint we impose an 

antisymmetry condition on the amplitudes of 

1 
1^) = ^ ^ijk^jkl), (1) 

ijkl=0 

as 

'^ijki = -'^kuj, (2) 

i.e. we impose antisymmetry in the first and second pairs of indices. 

An alternative (and more physical) way is the one of imposing such constraints on the 
original Hilbert space if ~ C^^ which renders to have a tensor product structure on one of 
its six dimensional subspaces 7i of the form 

H = (C'®C')A(C'®C2), (3) 

where A refers to antisymmetrization. As we know quantum tensor product structures are 
observable-inducedr^, hence in order to specify our system with a tensor product structure 
of Eq. ([3]) we have to specify the experimentally accessible interactions and measurements 
that account for the admissible operations we can perform on our system. For example we 
can realize our system as a pair of fermions with four single particle states where a part of 
the admissible operations are local unitary transformations of the form 

\^)^{U®V)^{U^V)\^), U,VeU{2), \^)eH. (4) 

Taken together with Eq. ([2]) this transformation rule clearly indicates that the first and 
the second and the third and fourth subsystems form two indistinguishable subsystems of 
fermionic type. 

The aim of the present paper is to study the interesting structure of the reduced density 
matrices of the form Qij, i < j, i,j = 1,2,3,4 arising from fermionic states that are 
elements of the tensor product structure as shown by Eq. ([3]). We can alternatively coin the 
term that these density matrices are ones with fermionic purifications. 

The organization of the paper is as follows. In Section II. we present our parametrized 
family of density matrices we wish to study. Using suitable local unitary transformations 



we transform this family to a canonical form. In Section III. based on these results we 
calculate the Wootters concurrence the negativity and the purity. We give a formula for the 
upper bound of negativity for a given concurrence. (We prove it in Appendix A.) In section 
IV. we analyze the structure of these quantities and discuss how they are related to each 
other. In particular we prove that the relevant entanglement measures associated with our 
four-qubit state satisfy the generalized Coffman-Kundu- Wootters inequality of distributed 
entanglemenlr^. For the residual tangle we derive an explicit formula, containing two from 
the four algebraically independent four-qubit invariants. In Section V. we investigate the 
Bures geometry of this special subclass of two-qubit density matrices. We show that thanks 
to our purifications being fermionic an explicit formula for the Bures metric with hyperbolic 
structure can be obtained. The conclusions and some comments are left for Section VI. 

II. THE DENSITY MATRIX 

Let us parametrize the 6 amplitudes of our normalized four qubit state |\E') of Eq. ([1]) 
with the antisymmetry property of Eq. ([2]) as 

"^ijkl = - i^ik^jl + ^ik^jl) , (5) 

where A and B are symmetric matrices of the form 

(zi - iZ2 -Zs \ (wi- iW2 -ws \ 

A= \ =e{zcr), B= \ =e{wcr), (6) 

y -z-i -zi - 1Z2J y -W3 -wi - 1W2J 

where z, w G C^, wcr = Wiai + W20"2 + ^3(73, with the usual ai Pauli matrices, 



(7) 



and the overline refers to complex conjugation. 

It is straightforward to check, that the normalization condition of the state |\1/) takes the 
form: 

ll^f = ||w||2 + ||z||^ := 1 (8) 

The density matrices we wish to study are arising as reduced ones of the form 

f? = f?i2 = Tr34(|^)(*|). (9) 




Notice that since the (12) and (34) subsystems are by definition indistinguishable we also 
have g = gi2 = Qsa- 

A calculation of the trace yields the following explicit form for g 

g = lil + A), (10) 

where 1 denotes the 4x4 identity matrix, and A is the traceless matrix 

A := xcr ® / + / ® ycr + wcr (g) zcr + wcr zcr, (11) 

x := zw X w, y := iz x z, (12) 

where / is the 2x2 identity matrix. Notice, that x, y G R^, and xw = xw = yz = yz = 0. 
Due to this, and the identities 

I |2 II ||4 2 — 2 I |2 II ||4 2—2 fio\ 

X = w — WW, y = z — zz, (13) 

II nil " it/i nil ' \/ 

it can be checked, that A satisfies the identity 

A' = il-r^')l, (14) 

where 

ri=\w^-z'^\, (15) 

< r/ < 1, (16) 

Notice that the quantity rj is just the Schliemann-measure of entanglement for two-fermion 
systems with 4 single particle states^^i^. Indeed our density matrix g (with a somewhat 
different parametrization) can alternatively be obtained^S as a reduced one arising from such 
fermionic systems after a convenient global f/(4), and a local U{2) x f/(2) transformation of 
the form / ® o"2 . 

Now by employing suitable local unitary transformations we would like to obtain a canon- 
ical form for g. According to Eq. (jlj) the transformations operating on subsystems 12 or 
equivalently 34 are of the form f/ (g) V" G f/(2) x f/(2). 

As a first step let us consider the unitary transformation 

f/(u, a) := e't"- = cos (|) ^ + ^ sin (^|) ucr. (17) 



which is a spin-i representation of an SU{2) rotation around the axis u G R^, (|u| = 1) 
with an angle a. A special rotation from x to x' (x' 7^ — x) can be written as 



[/(u, Q;)^(xcr)f/(u, a) = x'cr, 



U{u, a) 



v/2^ 



XXX 



XX 



= (x2/ + (x^)(xV)) 



,xx 



u 



a = arccos 



(18) 
(19) 

(20) 



|x X x'l ' ■ XX 

Employing this, we can rotate the vector x to the direction of the coordinate axis z. In 
this case 



r := X 



f/x 



^2r(r + X3 
[/t(xo-)f/, = ms. 



rl + (xcr)o-3) , 



(21) 
(22) 

(23) 



Moreover, using Eq. (iT9l) it can be checked that due to the special form of x (see Eq. ([T2l 
the transformation above rotates the third component of w into zero 



UI{w(t)U^ = w'o". 



wx 



w = w 



r^ + xx' 



[x + x'l 



^ r+X3 ■^ 



(24) 
(25) 



A similar set of transformations can be applied to ycr 



K: 



V^iy(T)Vy = sas. 

Kt(zCT)yy = ZV, 



{si + (ycr)cr3) 



5^ + yy' 



Zl 



s+ys 



Z3 



2 s+ys 3 



(26) 
(27) 

(28) 
(29) 

(30) 



Obviously, every U G U(2) unitary transformation acting on an arbitrary a G C'^ as 
Wacrll = a'cr preserves a^ and ||a|p, since a^ = — det(acr), and ||a|p = ^ Tr((acr)'''(acr)). 



Hence 



/2 2 

W = W , 



W 



/||2 



W 



/2 2 

z = z , 



,'l|2 



and 



7] = 1] 



(31) 
(32) 

(33) 



are invariant under local U{2) x f/(2) transformations. (The entanglement measure r] is also 
invariant under the larger group of U{4) transformations.) 

Now by employing the local U{2) x f/(2) transformations U^ Vy, our density matrix 
can be cast to the form, 



g' = (U^ ® \/y)t g (u^ ^ v;) = - (1 + A') , 
A' = ra^ ® I + I ® 5(73 + w'cr (g) z'cr + w'cr ® z'cr, 



(34) 
(35) 



where A' has the special form 



A' 



^3 ai — ia2 

/33 /9i-^/92 

Pi+il32 -Ps 

tti + ia2 —as 



(36) 



with the quantities ol^(3 E R^ defined as 



q: 



6-6 


C1 + C2 


r + s 



/3 



6 + 6 


Ci- 


-C2 


r - 


- s 



^2 = W'^l + ^i^24> 



Cl = ^2^1 + ^i^2^ 



1) 



C2 = W^Z^ + Wi^2- 



(37) 

(38) 
(39) 



Thanks to the special shape of A, we can regard g' as the direct sum of two 2x2 blocks, i.e. 
(/ + acr) and (/ + /3cr). Having obtained the canonical form of our reduced density matrix 
g, now we turn to the calculation of the corresponding entanglement measures. 



III. MEASURES OF ENTANGLEMENT FOR THE DENSITY MATRIX 

A. Concurrence 

In this section we calculate the Wootters-concurrence^ of our density matrix q defined in 
Eqs. f fTOj) - flT2l) . This quantity is defined as 

C = max{0, Ai - A2 - A3 - A4} (40) 

where Ai > A2 > A3 > A4 are the square roots of the eigenvalues of the matrix- gg where 



g = {a2 ® (J2) g{(J2 ® (J2) . (41) 

This matrix (the Wootters spin-fiip of g) is known to have real nonnegative eigenvalues. 
Moreover, the important point is that C is an SL{2, C) x SL(2, C) invariant^, hence we 
can use the canonical form we obtained in the previous section via using the transformation 
f/x ® 1/y G SU{2) X SU{2) C SL{2, C) x SL{2, C) for its calculation. 

It is straightforward to check that 16g'g' has the same X-shape as g' , with the blocks 
{aol + acr) and {PqI + /3cr) where a^, [3y n^u = 0, 1, 2, 3 are quadratic in a. j3: 



a, 



1 + al + ai 



an 



P. 



2a\ — i2a2aj, 

2a2 + i2an,a\ 



1 + /3? + /3| - /3| 

2/?i - 22/52/?3 
2/?2 + ^2/?3/?l 





gC^ 



(42) 



gC^ 



^2 



The eigenvalues of the blocks {(XqI + acr) and {J5qI + ^cr) are a^ ± va^ and /3o ± \//3' 
respectively. Now, we can express these with the help of the quantities a, /3 of (l42l) and get 
the eigenvalues of (J (J in the form 



A^, 




a\\ a\^\j\ — 



(43) 



Now, using Eqs. (1371) - (139|) . we have to express these as functions of our original quantities 
z^, w^, ||z|p and ||w|p. Straightforward calculation shows, that: 

al + al = 2||w'f llz'f + w'V^ +^'2^/2 _ 3^^ 

j3l + (3l = 2||w'f llz'f + w'V2 + w'2z'2 + 2rs, 

(44) 

1 2 oil /||2l|/||2i /2 — /2 I /2— /2 o 

1 — ctg = z||w ||||z ll+ww+zz — 2rs, 

1 /92 oil / II 2 II /||2 , /2 — /2 I /2— /2 , o 

1 — P3 = 2||w II ||z II + w w + z z + 2rs. 

The formulas above are expressed in terms of quantities invariant under our transformation 
yielding the canonical form (see Eqs. fl3T|) - fl32|) ). hence we can simply omit the primes. Hence 
by using Eq. (iT3l) and (ITSll we can establish that 



(45) 



2 I 2 1 2 2 

For further use, denote: 

7_^ := r + s = ag, 7. := r - s = P3. (46) 

With these, the square root of the eigenvalues of gg are 




1 



5 



4 , v1-7-±a/1-7--T 



(47) 



The biggest one of these is Xmax = \ ( \^^ — 7- + a/1 — 7^ — r^^ j and after subtracting the 
others from it, we get finally the nice formula for the concurrence 



C{g) = max |o, ^ (^1 - 7' - ^^ - \A^) } (48) 

with the quantities defined in Eqs. ( TT2l) . ( TTSj) . ( 12T1) . f l26l) and ( H6|) containing our basic 
parameters w and z of f). One can easily check by the definitions fH6l) . that the surface 
dividing the entangled and separable states in the space of these density matrices is a special 
deformation of the rj = Klein-quadric,— given by the equation: 

77^ = 4rs. (49) 

This can be also seen from the (l52l) formula of negativity, see in next subsection. 



B. Negativity 

Another entanglement-measure which we can calculate for q is the negativity. It is related 
to the notion of partial transpose and the criterion of Peres^^. It is defined by the smallest 
eigenvalue of the partially transposed density matrix, as foUows^*^ 



J\f{Q) = max {0, -2fLmin} ■ 



(50) 



Since the eigenvalues of a complex 4x4 matrix are invariant under U{A) transformation, 
we can use again the SU{2) x SU{2) C f/(4)-transformed g' of Eq. (!M|) . 

Denote by g''^'^ the partial transpose of g' with respect to the second subsystem. This 
operation results in the transformation A'^^ = rcr3C>?>/ + /(8>S(T3 + w'<T® (z'cr)-^ + w'<T® (z'cr)-^ 
i.e. only z'2 changes to —z'2. By virtue of this, retaining the (l38l) . fl39l) definitions of ^1, ^2, 
6, C25 and redefining a, /3 G R'^ of Eq. fl37|) as 



a 



6 + 6 




6-6 


C1-C2 


, /3 = 


C1 + C2 


r + s 




r — s 



(51) 



the calculation proceeds as in section [Til A[ The eigenvalues of g^"^ are /Xj = {|(l±|a|),i(l± 
|/3|)}, and straightforward calculation shows, that ol^ = 1 — rf + Ars and (3 = 1 — rf — Ars. 
Hence one can see, that the negativity of g is 

M{g) = max |o, ^ (^1 - 7^2 + 4rs - 1 U , (52) 

with the usual conventions of equations (IT5|) . fl2Tl) and fl26j) . 



C. Comparsion of concurrence and negativity 

For a 2-qubit density matrix we can write the following inequalities between the concur- 
rence and the negativity 



v/(l-C)2 + C2 - (1 - C) < AT < C, 



(53) 



which are known from a paper of Audenaert et. al.^ Our special case with fermionic corre- 
lations may give extra restrictions between concurrence and negativity, so we can pose the 
question, whether we can replace inequality (l53l) by a stronger one. 

10 




FIG. 1: Range of value of negativity for a given concurrence, with its boundaries, as in (jMj) . 



Fig. [I] shows the resuh of a numerical calculation. The gray field denotes the possible 
entangled states of the (TTU1) -form on the C — Af plane. The upper bound of these can be 
analitically determined, and it can be proven, (see Appendix |A|) that the following inequality 
holds for A/": 



v/(l-C)2 + C2 - (1 - C) < AT < i (^V2-(1-2C)2 - l) < C 



(54) 



To see, that this upper bound is the tightest, consider the special case, when w = z. These 
states realize the boundary, so the second inequality in (IM|) turns to equality. (In this 

0, and for entangled states, C 



case 7] 



0, r 



s, 7+ 



2r, 7_ 



1(1 



Vl-4r2), 



A/" = ^(Vl + 4r2 — 1). These depend only on r, wich can be expressed from C, thus we can 
express the negativity of these states with their concurrence, and get back the curve of the 
upper bound.) 

It can be seen, by calculating the intersection of the corresponding curves of (IMll . that for 
maximally entangled states C{gmax) = |, and Mi^Qmax) = 2^ ■ We can study the behavior 
of these states: from the (HHI) concurrence formula one can see, that 



C 



[if = 0, and 7! = 0, and 7^ = l) 



(55) 



The first two of these constraints hold, if and only if w^ = z^, and r = s, because of ( TT5l) . 



T3l) and fl46l). If w^ = z^ then r = s and llwlP 



I are equivalent, and if r = s then 



7^ 



4^2 



\w\ 



W 



2|2 



T and follows, that w^ = 0. 



C = C 



1 
2 



w 



-, and w^ 



(56) 



11 



Since the ( !M|) transformation on q preserves the quantities appearing here, we can easily 
calculate the flM|) canonical form of the maximally entangled state q' . Let us choose an ansatz 



iT 



and 



of the form ^ and ^ for w'^^^ and z'„„^ as w'^^^ = -^ 

r 1^ 

'^'max = 73 cos(/5)e*'^^,sin(/3)e*'^^0 • These satisfy the first constraint of ([56]), and from 
the second follows that cos a = sina = 4= and (pi = {p2 — ^ =: (p and the same for z'^^^x- 



w^ 



a*^ 



3# 



(57) 



Then for the density matrix with maximal concurrence we get the expression 



1 
4 



§00 


o" 


1 e'^ 





e-'^ 1 








1 

2. 



(58) 



with 6 = (f—ip is the only parameter characterizing this maximally entangled density matrix. 



D. Purity 



The purity is measuring the degree of mixedness of a density matrix. For our g thanks 
to the special property of A (see in Eq. flT^ ) it can easily be calculated. We have the result 

4 - ^ - 2 
by virtue of Eq. flT6l) . The participation ratio is by definition 

1 4 



R 



Tr ^^ 2 — T]'^ 



(59) 
(60) 

(61) 



IV. RELATING DIFFERENT MEASURES OF ENTANGLEMENT 



Now we would like to discuss the physical meaning of our quantities derived in the previous 
section. First of all let us notice that the 



^1 = Tr234(|^)(^|) = TT2{gi2) = TT2{g) = ^{I + ^(t) 



(62) 



12 



Q2 = Tri34(|^)(^|) = Tri(f?i2) = Tri(f?) = ^(/ + ycr) (63) 

reduced density matrices describe the entanglement properties of subsystems 1 and 2 to the 
rest of the system described by the four-qubit state |\1/). It is well-known that the measures 
describing how much these subsystems are entangled with the rest are C^^m) ~ 4Det(^i) 
and ^2^34) = 4Det(^2)- Due to Eqs. fl62tl63p these quantities are 

< Ci%34) = 1 - r^ < 1, < C|(i34) = l-s^<l. (64) 

Clearly, since ^i = ^3 and §2 = Qi, we also have 

^3(124) = ^1(234) 5 ^4(123) ~ ^2(134)- ("5) 

Moreover, we already know that gi2 = ^34 = ^. A straightforward calculation of the 
two-partite density matrices f)i4 and f)23 shows that they again have the form of Eq. (10) 
with the sign of w is changed in the first case and the vectors w and z are exchanged in the 
second. Since these transformations do not change the value of the concurrence, we have 

^12 — ^14 — ^23 ~ ^34- (.66) 

Now the only two-qubit density matrices we have not discussed yet are the ones ^13 and ^24- 
Their form is 

(^13)iA:j'fc' = 7^ (ll^^ll ^ik^i'k' + ISikBi'k') , i024)jij'i' = " (||w|| Sjl^j'l' + AjlAjn') . (67) 

Recall now the that the (jlj) transformation property of our four-qubit state gives rise to the 
corresponding ones for the reduced density matrices 

gi3^{U^U)gi3(U®U), §24^ {V ^V)g2i(V ^V). (68) 

For U,V E SU{2) we have VeV* = UeU* = e, hence the tensors occurring in Eq. ( l67|l 
transform as 

e^e, A^VAV\ B^UBU'. (69) 

Using the (j6]) definition of A we have for example 



VAV = eVzWV' = eVzaV^ = ez'a- = ez'W, (70) 



13 



where by choosing V = V^ oi Eq. (27) we get for z' the (30) form. Finally these manipula- 
tions yield for ^24 the canonical form 

' Kq + K3 Hi — iK2^ 





^24 



IwlP -llwlP 



-llwlP llwlP 



\Ki + iK2 



Kq- K3J 



where 



Ko = 

Notice that 



'l|2 



-Ki 



I ' |2 

\z, — \Z. 



J\2 
2I 7 



Kn 



nji ~t~ njo ~t~ njn 



(71) 



-K2 = 2Re(zi4), -/^s = 2Im(2i4). (72) 



(73) 



hence the eigenvalues of ^24 are ||w|p, ||z|p,0,0, i.e. our mixed state is of rank two. The 
structure of qi^ is similar with the roles of w and z exchanged. Following the same steps as in 
Section III. A. we get for the corresponding squared concurrences the following expressions 

2|n2 ^9 /ii ii9 i9i\2 



'-'IS 



w 



'\r 



'-'24 



W 



1) • 



(74) 



Let us now understand the meaning of the invariant r] from the four-qubit point of view. 
It is known that we have four algebraically independent SL{2, C)®^ invariants^ii^ denoted 
by H, L, M and D. These are quadratic, quartic, quartic and sextic invariants of the complex 
amplitudes "^ijki respectively. The invariants H and L are given by the expressions 



H = ^0^15 - ^1^14 - ^2^13 + ^3^12 - ^4^11 + ^5^10 + ^6^9 " ^7^ 



7^8, 



(75) 



and 



Det 



(76) 



\I'4 \l/5 \I'6 "^7 
\*12 ^13 ^14 ^15/ 

where instead of the binary one we used the decimal labelling. For the explicit form of the 
remaining two invariants M and D see the paper of Luque and Thibon^^. A straightforward 
calculation shows that for our four-qubit state we have M = D = however, 



H = -^{z^ + ^') 



1 

16" 



:z2-wY, 



(77) 



14 



hence 



1^1 = ^V'. (78) 



For convenience we also introduce the quantity 

a= Iw^ + z^l =2|i7|. (79) 

Hence t] = |w^ — z^| and a = |w^ + z^| are related to the only nonvanishing four qubit 
invariants L and H. Using the definitions of these quantities and Eq. (13) one can check 
that 

^•2 2 I / 2 I 2\ oil 2||| 2| ^2 2 , / 2 , 2\ oil 2||| 2| ron\ 

W3 = ■s + 0(7 +cr ) ~ 2||z^|||w^|, C24 = r^ + -(77 +(T ) - 2||w^|||z^|. (80) 

Hence we have the inequality 

C!, + Cl<s' + r' + v' + ^'. (81) 

Moreover, since 0^2 = C\^ after taking the square of Eq.(48) we get 



CI2 + Cl = l-r^-s'--r^'- ^(1 - r^2 _ ^2 )(i _ ^2 ) (53) 

Combining this result with Eqs.(64) and (80) we obtain 

^12 + ^3 + ^14 + ^1 = W(234) C^2 + ^23 + ^24 + ^2 = ^3(134), (83) 



where 



Si = 2||z||V'l + \l {\^' + p+){\^' + P-) - \y\ (84) 



Here 



S2 = 2||w|nz2| + ^/(^or2 +p+)(la2 +p_) - \a\ (85) 



p± = 2||z|n|w|p±-(4rs-r7^). (86) 



Notice that by virtue of Eq.(13) p_ is nonnegative. Moreover according to Eq. (52), for 
nonseparable states (^12, ^14, f?34, P23) we have nonzero negativity hence 4rs > rf hence p+ 
is also nonnegative. In this case the residual tangles Si and S2 as defined by Eqs. (84-85) 
are positive as they should be, hence the generalized Coffman-Kundu-Wootters inequalities 
of distributed entanglement ^^'^* hold 

^12 + ^13 + ^14 — ^1(234) ^12 + ^3 + ^24 — ^(134)- (87) 
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For separable states we have Cu = Cu = Cu = C23 = and a calculation shows that the 
(87) inequalities in the form Cf^ < C^(234) ^^^ ^24 — ^2(134) ^^^^^ ^^^"^ with residual tangles 

v-i oil 112/1 2| I II ||2\ v-i oil 112/1 2| , || ||2\ /oo\ 

Li = 2||z|| (|w I + ||w|| ) L2 = ^||w|| (|z I + ||z|| ). (88) 

Eqs. (84), (85) and (88) show the structure of the residual tangle. Unlike in the well- 
known three-qubit case these quantities among others contain two invariants 77 and a char- 
acterizing four-partite correlations. The role of a (which for a general four-qubit state is a 
permutation-invariant) is to be compared with the similar role the permutation invariant 
three-tangle T123 = 4|P| plays (an SL{2, C)®^ invariant) within the three-qubit context. {T> 
is Cayley's hyperdeterminanti^.) An important difference to the three-qubit case is that the 
residual tangles Si 2 seeem to be lacking the important entanglement monotone property. 
However, according to a conjecture the sum Si + S2 could be an entanglement monotone. 
We hope that our explicit form will help to settle this issue at least for our special four-qubit 
state of Eqs. (1-2). 

V. BURES METRIC 

As we have emphasized our density matrix g can be regarded as a reduced density matrix 
of a two-particle system on (C^ ® C^) A (C^ ® C^), meaning 







'^¥, (89) 



where \1/ is the 4x4 antisymmetric matrix occurring in Eq. ( 1761) . In the space of such 
fermionic purifications of our density matrix the curve \i/(t) is horizontal, when the differen- 
tial equation 

^t^ = ^t^ (90) 

holds. We can satisfy this equation by the ansatz 

^ = G^, G = G^ (91) 

for some G = G\ so that 

dg = {G, g}. (92) 

We can define the Bures metric on the space of density matrices, as follows^ 

dsl = -TT{Gdg). (93) 
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Let us now take into account the condition \1/-^ = — \1/. Taking the transpose of Eq. (1901) . 

we get 

(94) 

Using this result, we get a simpler formula for dg: 

dg = ^^t + ^^t = 2G^^^ = 2Gg. (95) 

and to get G, we only have to invert g. A calculation of the eigenvalues of g shows that they 
are of the form^ 

A± = 1(1 ± x/T^). (96) 

Hence g is nonsingular iff 77 7^ (i.e. iff |L| 7^ 0). In the following we consider this case. 
For nonsingular density matrices by virtue of Eq. flT4|) . g^^ can be calculated easily 

g-' = ^,{l-A), (97) 



hence: 



and the Bures-metric: 



G = -dgg-^ = ^{dA- dAA) , (98) 



1 



dsl = T^ Tr {dAdA - dAAdA) . (99) 

Since A is idempotent and traceless, one can see, that the trace of the second term equals to 
zero: 2TT{dAAdA) = Tr((iA(i(A^)) = Tr((iA(i(— r/^)!) = 0. Let us introduce the quantities 
fij = w.jZj + WiZj. With this notation we have 

ds% = — {dxidxi + dyjdyj + dfijdfij) , (100) 

(summation on i,j = 1, 2, 3 is implied.) Moreover, a calculation shows that 77^ = 1 — (xjXj + 
Hji/j + fijfij), so after putting the quantities Xj, yj, fij into a 15 component vector k G R^^ 
our final result is the nice formula 

4rfk2 



d^^ - ^ dk^ - -r? 
ass - ;l -k2 ~ A 
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k2)^ 



(101) 



Let us compare this formula with the one obtained for the Bures metric of one-qubit 
density matrices arising as a reduced density matrix from a pair of distinguishable qubits^S 

dll = ^^ (dp^ + sinh^ pdn^) , (102) 

4 cosh /? 

where 1 — C^ = tanh /? with C is the pure state concurrence for two qubits, and dn"^ is 
the usual line element on the two-sphere 5*^ expressed by the angular coordinates i) and 
(fi. Since the space of one-qubit density matrices is the Bloch-ball B'^ this parametrization 
provides a map between the upper sheet of the double sheeted hyperboloid H^ and B^. The 
standard metric on H^ is just dp"^ + smh"^ pdil"^. Hence we see that the Bures metric is up 
to the conformal factor C^/A is just the standard metric on the upper sheet of the double 
sheeted hyperboloid H^. However, using the stereographic projection one can show that 

dp' + sinh^ pdn' = 7--^, (103) 

(1 — K,"^)^ 

where -R1-R2 and R3 can alternatively be used to parametrize B^. Hence we can write 

4dR2 



dll = -C 

4 



1-R 



2^2 



(104) 



where the metric on the right is the standard Poincare metric on the unit ball which is now 
just the Bloch-ball. Comparing this equation with our previous expression of Eq. (1 10 11) we 
see that it is up to the conformal factor rj'^/A is just the Poincare metric on the Poincare 
ball B^^. We emphasize however, that unlike the usual one-qubit mixed state where all the 
Bloch parameters characterizing the density matrix are independent, here the 15 parameters 
associated to the vector k are subject to nontrivial constraints. These constraints describe 
some nontrivial embedding of the space of nonsingular density matrices V into the Bloch 
ball B^^ with our Bures metric of Eq. (101). 

VI. CONCLUSIONS 

In this paper we investigated the structure of a 12 parameter family of two-qubit den- 
sity matrices with fermionic purifications. Our starting point was a four-qubit state with 
a special antisymmetry constraint imposed on its amplitudes. Such states are elements 
of the space (C^ ® C^) A (C^ C^) and the admissible local operations are of the form 
{U ® V) <S) {U <S) V) G SL{2, C)®'*. Our density matrices are arising as the reduced ones 
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Q = Qi2 = Tr34(|^)(\l/|). Since the 12 subsystem is indistinguishable from the 34 one we 
have gi2 = Qsa- We obtained an exphcit form for g in terms of the 6 independent complex 
amplitudes w and z of our four-qubit states. Employing local unitary transformations of 
the form U ®V e SU{2) x SU{2) C SL{2, C) x SL{2, C) we derived the canonical form 
for g. This form enabled an explicit calculation for different entanglement measures. We 
have calculated the Wootters concurrence, the negativity, and the purity. The quantities 
occurring in these formulae (and some additional ones) are subject to monogamy relations 
of distributed entanglement similar to the ones showing up in the Coffman-Kundu- Wootters 
relations for three-qubits. They are characterizing the entanglement trade off between dif- 
ferent subsystems. We have entanglement measures rj and a describing the intrinsically 
four-partite correlations, quantities ( Wootters concurrences) keeping track the mixed state 
entanglement of the bipartite subsystems embedded in the four-qubit one. Finally we have 
the independent quantities ^^234-) and ^2^34-) measuring how much subsystems 1 and 2 are 
entangled individually to the rest. We derived explicit formulas displaying how these im- 
portant quantities are related. At last we have studied in some detail the Bures geometry 
underlying the structure of these density matrices. We have shown that the constraint 
of antisymmetry makes it possible to obtain a nice explicit formula for the Bures metric 
reminiscent of the ones known from hyperbolic geometry^. 
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APPENDIX A: UPPER BOUND OF NEGATIVITY 

In this fermionic-correlated case, defined by equations (ITU|) . (ITT]) . (IT^ and (ITB]) . we can 
prove the following inequality: 
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Theorem: For all entangled q: 

U{q) < \ (v/2-(l-2C(^))2 - l) . (Al) 

Proof: Insert Eqs. (gH]) and ([52]) into flXTjl : 



1 . . , , . , , x2 



2 vvl-^' + 4rs-lj<2 I V2-(l-\/l-7'-^' + \/l-7i) -1| , (A2) 



after some algebra, we can rearrange the terms: 



< 2r/^ - 2 + 7^ + 7^ - 4rs 



+ 2^1-7^-r/2-2^1-7i + 2^1-7^-r/yi-7i 
If follows from the definition fH6|) that 7^ — 4rs = 7^ . With this: 



(A3) 



The right-hand side is factorizable: 



0< (^^l-7^-r/2-^l-7iJ \\-^j\-^l-ri^y (A5) 

The second parenthesis is obviously positive. For entangled states C{^q) > 0, and the first 
parenthesis is proportional to the concurrence, which is strictly positive. Q.E.D. 
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